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Abstract. We give a short review of hydrodynamic models at heavy ion collisions 
F ^ , from the point of view of initial conditions, an equation of states (EoS) and freezeout 

[T i ' process. Then we show our latest results of a combined fully three-dimensional 

■ macroscopic/microscopic transport approach. In this model for the early, dense, 

£NJ \ deconfmed stage relativistic 3D- hydrodynamics of the reaction and a microscopic non- 

equilibrium model for the later hadronic stage where the equilibrium assumptions are 
not valid anymore are employed. Within this approach we study the dynamics of hot, 
' bulk QCD matter, which is being created in ultra-relativistic heavy ion collisions at 

OO RHIC. 
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1. Hydrodynamic Models at RHIC 

The first five years of RHIC operations at ^/snn = 130 GeV and ^/snn = 200 GeV 
have yielded a vast amount of interesting and sometimes surprising results. There exists 
mounting evidence that RHIC has created a hot and dense state of deconfined QCD 
matter with properties similar to that of an ideal fluid [fl [2] - this state of matter has 
been termed the strongly interacting Quark- Gluon- Plasma (sQGP). One of the evidence 
is the success of ideal hydrodynamic models in various physical observables. Especially 
for the first time, in elliptic flow the hydrodynamic limit shows good agreement with 
the experimental data at RHIC, though at AGS and SPS hydrodynamic models give 
larger value as compared to experimental data. 

The sophisticated 3D ideal hydrodynamic calculations, however, reveal that many 
of experimental data have not yet been fully evaluated or understood [5]. For 
example, elliptic flow at forward and backward rapidity [j] and at peripheral collisions 
is overestimated by ideal hydrodynamical models. The Hanbury-Brown Twiss (HBT) 
puzzle is not completely understood from the point of view of hydrodynamics. These 
distinctions between hydrodynamic calculation and experimental data suggest that the 
ideal hydrodynamic picture is not applicable to all physical observables even in low 
transverse momentum region. At present our main interest in hydrodynamic models is 
the following: How perfect the sQGP is? 

| Brazil group shows improved results of V2 at forward/backward rapidity by using event-by-event 
fluctuated initial conditions [3]. 
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Here, first, we shall give a short review of hydrodynamic models at heavy ion 
collisions. In addition to a numerical procedure for solving the relativistic hydrodynamic 
equation, hydrodynamic models are characterized by initial conditions, EoSs and 
freezeout process. 

The necessity of input of initial conditions for hydrodynamic models is one of the 
largest limitations of them. Because an initial condition is not be able to determined 
in the framework of hydrodynamic model itself, usually a parametrization of energy 
density and baryon number density based on Glauber type is used and parameters in 
it are determined by comparison with experimental data [H [11]. Recently there are 
some studies in which more basic approaches, Color Glass Condensate (CGC) (6], pQCD 
+ saturation model [7] are used for construction of initial conditions. 

The most important advantage of hydrodynamic models is that it directly 
incorporates an EoS as input and thus is so far the only dynamical model in which a 
phase transition can explicitly be incorporated. In the ideal fluid approximation - and 
once an initial condition has been specified - the EoS is the only input to the equations 
of motion and relates directly to properties of the matter under consideration. In this 
sense a hydrodynamic model is a bridge between QCD theory and experimental data and 
indispensable to describe heavy ion physics. However in usual practical hydrodynamic 
simulations, an EoS with 1st order phase transition (Bag model) is used. In fact, there 
are few studies on effect of order of QCD phase transition on physical observables [8J. 

Conventional hydrodynamic calculations need to assume a freezeout temperature 
at which the hydrodynamic evolution is terminated and a transition from the zero 
mean-free-path approximation of a hydrodynamic approach to the infinite mean-free- 
path of free streaming particles takes place. The freezeout temperature usually is a 
free parameter which can be fitted to measured hadron spectra. There are several 
approaches for dealing with freezeout process: chemical equilibrium [H [8], partial 
chemical equilibrium [5], continuous emission model [9] and construction of a hybrid 
model of a hydro + cascade model [TUJ [TlJ [12] . 

In Tab. [1] several hydrodynamic models are listed from the viewpoint of initial 
conditions and freezeout processes, because almost the same EoS with strong 1st order 
phase transition is used. Reference [4] presents the first calculation which shows the 
remarkable agreement with experimental data for both of Pt spectra and V2 at RHIC. 
However it turns out that the assumption of single freezeout temperature where chemical 
freezeout and kinetic freezeout occur at the same time fails in reproducing hadron ratios 
correctly. To obtain correct proton Pt spectra, we need to renormalize the Pt spectra 
using the p to n ratio at the critical temperature. In addition, recently, Hirano and 
Gyullasy point out that the good agreement of elliptic flow with experimental data 
may accidentally happen. Hirano and Kolb et al. propose that the introduction of two 
kinds of freezeout processes, chemical freezeout and kinetic freezeout to hydrodynamic 
models. In this model, normalization of the Pt spectra for each particle are obtained 
correctly [5], but to get better agreement with experiments in elliptic flow additional 
initial transverse flow is needed. At present, the combination of Gluaber type and a 
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cascade model gives us the most promising result for both of Pt spectra and v-i- Hirano 
et al. perform calculations using CGC for initial condition and a cascade model for 
freezeout process, which suggests that viscosity is not negligible even at early stage 
of the expansion. In the next section we show our latest results [11] based on the 
full 3D hydrodynamic approach [13] with the microscopic Ultra-relativistic Quantum- 
Molecular-Dynamics (UrQMD) model [14J. 



Table 1. Hydrodynamic models at RHIC. 



Reference 


Initial Conditions 


Freezeout Process 


i 


Glauber type 


chemical equilibrium 


m 


Glauber type 


partial chemical equilibrium 




Glauber type 


cascade model 


m 


CGC 


partial chemical equilibrium 


LX2J 


CGC 


cascade model 



2. Hydro+UrQMD Model 



We calculate hadron distribution at switching temperature from the 3D hydrodynamic 
model using Cooper- Frye formula [16] and produce initial conditions for UrQMD model 
by Monte Carlo from it. Such hybrid macro/micro transport calculations are to date 
the most successful approaches for describing the soft physics at RHIC. The biggest 
advantage of the hydrodynamic description is that it directly incorporates an EoS as 
input - one of its largest limitations is that it requires thermalized initial conditions and 
one is not able to do an ab-initio calculation. 

Figure [H shows a schematic sketch of the full 3D hydrodynamic model + UrQMD. 
After heavy ion collisions, first, hydrodynamic expansion starts. We introduce one 
more parameter, switching temperature from hydrodynamic picture to hadron base 
event generator, UrQMD. This switching temperature should be just below the critical 
temperature. Here it is set to 160 MeV. 



Full 3D Hydrodynamics 



hadronization 



EoS:1st order phae transition 4 9° ope f~ Frye 
r : formula 




UrQMD 



final state 
j Monte Carlo interactions 

7c 7sw ffm 



Figure 1. Schematic sketch of 3D hydro+UrQMD model. T c (= 160 MeV) and T sw (= 
150) MeV are critical temperature and switching temperature from hydrodynamics to 
UrQMD model, respectively. 
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Figure 2 shows the Pt spectra of 7r + , K and p at a/s/vw = 200 GeV central 
collisions. The most compelling feature is that the hydro+micro approach is capable of 
accounting for the proper normalization of the spectra for all hadron species without any 
additional correction as is performed in the pure hydrodynamic model. The introduction 
of a realistic freezeout process provides therefore a natural solution to the problem of 
separating chemical and kinetic freeze-out in a pure hydrodynamic approach. 

In Fig. 3 centrality dependence of Pt spectra of 7r + is shown. The impact parameter 
for each centrality is determined simply by the collision geometry. The separation 
between model results and experiment appears at lower transverse momentum in 
peripheral collisions compared to central collisions, just as in the pure hydrodynamic 
calculation. The 3D hydro + micro model does not provide any improvement for this 
behavior, since the hard physics high Pt contribution to the spectra occurs at early 
reaction times before the system has reached the QGP phase and is therefore neither 
included in the pure 3D hydrodynamic calculation nor in the hydro+micro approach. 
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Figure 2. Pt spectra for tt + , 
K + and p at central collisions with 
PHENIX data Q7|. 



Figure 3. Centrality dependence 
of P T spectra of vr+ with PHENIX 
data Q2]. The P T spectra at 10-15 
%, 15-20 % and 20-30 % are divided 
by 5, 25 and 200, respectively. 



Figure 4 shows the centrality dependence of the pseudorapidity distribution of 
charged hadrons compared to PHOBOS data [18]. Solid circles stand for model results 
and open circles denote data taken by the PHOBOS collaboration [18]. The impact 
parameters are set to b = 2.4,4.5,6.3,7.9 fm for 0-6 %, 6-15 %, 15-25 % and 25- 
35 % centralities, respectively. Our results are consistent with experimental data 
over a wide pseudorapidity region. There is no distinct difference between 3-D ideal 
hydrodynamic model and the hydro + UrQMD model in the centrality dependence of 
the psuedorapidity distribution, indicating that the shape of psuedorapidity distribution 
is insensitive to the detailed microscopic reaction dynamics of the hadronic final state 
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In Fig. 5 we analyze the Pt spectra of multistrange particles. Our results show good 
agreement with experimental data for A, S, Q for centralities 0-5 %. In this calculation 
the additional procedure for normalization is not needed (Fig. 6). Recent experimental 
results suggest that at thermal freezeout multistrange baryons exhibit less transverse 
flow and a higher temperature closer to the chemical freezeout temperature compared 
to non- or single-strange baryons [191 [20] . This behavior can be understood in terms 
of the flavor dependence of the hadronic cross section, which decreases with increasing 
strangeness content of the hadron. The reduced cross section of multi-strange baryons 
leads to a decoupling from the hadronic medium at an earlier stage of the reaction, 
allowing them to provide information on the properties of the hadronizing QGP less 
distorted by hadronic final state interactions 



Charged Hadrons 
O PHOBOS (0-6 %, 6-15 %, 15-25 %, 25-35 %) 
• Hydro + UrOMD (b=2.4, 4.5, 6.3, 7.9 fm) 
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Figure 4. Centrality dependence 
of pseudorapidity distribution of 
charged particles with PHOBOS 
data [TS]. 



Figure 5. Pt spectra of multi- 
strange particles at centralities 0- 
5 % and 10-20 % with STAR data 



In Fig. 7 the mean transverse momentum (Pt) as a function of hadron mass is 
shown. Open symbols denote the value at T sw = 150 MeV, corrected for hadronic decays. 
Not surprisingly, in this case the (Pt) follow a straight line, suggesting a hydrodynamic 
expansion. However if hadronic rescattering is taken into account (solid circles) the 
(Pt) do not follow the straight line any more: the (Pt) of pions is actually reduced 
by hadronic rescattering (they act as a heat-bath in the collective expansion), whereas 
protons actually pick up additional transverse momentum in the hadronic phase. RHIC 
data by the STAR collaboration is shown via the solid triangles - overall the proper 
treatment of hadronic final state interactions significantly improves the agreement of 
the model calculation with the data. 

In Fig. 8 we plot elliptic flow v<i as a function of Pt- The solid line stands for the 
pure hydro calculation, terminated at the switching temperature T sw and solid circles 
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Figure 6. Pt spectra for multi- 
strange baryons at central collisions 
from pure hydro with STAR data 
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Figure 7. Mean Pt as a function 
of mass with STAR data (Au+Au 
^sWn = 130 GeV) [20]. 



denote the full hydro+micro calculation. We find that the QGP contribution to the 
elliptic flow depends on the transverse momentum - for low nearly 100% of the 
elliptic flow is created in the QGP phase of the reaction, whereas the hadronic phase 
contribution increases to 25% at a Pt of 1 GeV/c. 

Figure 9 shows the elliptic flow as a function of r/: the pure hydrodynamic 
calculation is shown by the solid curve, the hydrodynamic contribution at T sw is denoted 
by the dashed line and the full hydro+micro calculation is given by the solid circles, 
together with PHOBOS data (solid triangles). The shape of the elliptic flow in the pure 
hydrodynamic calculation at T sw is quite different from that of the full hydrodynamic 
one terminated at a freeze-out temperature of 110 MeV. Apparently the slight bump at 
forward and backward rapidities observed in the full hydrodynamic calculation develops 
first in the later hadronic phase, since it is not observed in the calculation terminated 
at T sw . Evolving the hadronic phase in the hydro+micro approach will increase the 
elliptic flow at central rapidities, but not in the projectile and target rapidity domains. 
As a result, the elliptic flow calculation in the hydro+micro approach is closer to the 
experimental data when compared to the pure hydrodynamic calculation. 

We developed a novel implementation of the well known hybrid macro- 
scopic/microscopic transport approach, combining a newly developed relativistic 3+1 
dimensional hydrodynamic model for the early deconfined stage of the reaction and the 
hadronization process with a microscopic non-equilibrium model for the later hadronic 
stage. 

Within this approach we have dynamically calculated the freezeout of the hadronic 
system, accounting for the collective flow on the hadronization hypersurface generated 
by the QGP expansion. We have compared the results of our hybrid model and of 
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Figure 8. Elliptic flow as a 
function of Pt of 7r + at centrality 
5-10 % with STAR data [21]. 



Figure 9. Elliptic flow as a 
function of rj of charged particles 
with PHOBOS data [22]. 



a calculation utilizing our hydrodynamic model for the full evolution of the reaction 
to experimental data. This comparison has allowed us to quantify the strength of 
dissipative effects prevalent in the later hadronic phase of the reaction, which cannot be 
properly treated in the framework of ideal hydrodynamics. 

3. Summary 

The full 3D relativistic hydrodynamics + cascade model is one of successful and realistic 
models for description of dynamics of hot QCD bulk matter at RHIC, which helps us to 
understand medium property at RHIC in detail. Using this model, we can explore 
interesting phenomena which are caused by interactions between medium and jets 
[23] EI]. One of proposed interesting physical observables is mach cone [24] for the 
wake of them, from which we can know medium property in detail, too. 

However still further investigation for EoS and initial conditions in hydrodynamics 
is needed. Especially the EoS is the key to know the QCD phase transition directly from 
comparison with experimental data. Recent lattice QCD calculation shows relatively 
high critical temperature with crossover phase transition and the existence of critical 
end point [25], which suggests that the EoS with strong 1st phase transition is not 
realistic. Becuase in hydrodynamic calculations, however, outputs may also be changed 
easily by choice of initial conditions and freezeout process, which makes it difficult to 
obtain the conclusive discussion on the EoS. First of all, we have to determine the 
most realistic initial conditions and freezeout process in hydrodynamic models before a 
detailed discussion about the EoS. 

Furthermore, viscosity effect in medium also starts to be discussed actively |26j . 
Due to serious difficulty in construction not only of a viscous hydrodynamic code 
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but also of framework of viscous hydrodynamics without the causality problem, the 
progress of study of viscous hydrodynamics has been slow. However, by virtue of recent 
rapid development, the practical calculations with viscous hydrodynamics which are 
comparable to experimental data will be achieved in the near future |26j. 

Finally, at LHC heavy ion collisions experiments will start in a year. Upcoming 
LHC data shall bring a lot of interesting, fruitful and even unexpected results for QGP 
physics, where hydrodynamic description will be helpful and useful for understanding 
of it. 
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